Metallic atomic wires on a patterned dihydrogeneted Si(OOl) 
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Electronic structure calculations for atomic wire of metals like Al, Ga and In are performed 
for a patterned dihydrogeneted Si(001):l x 1 in search of structures with metallic behavior. The 
dihydrogeneted Si(001) is patterned by depassivating hygrozen atoms only from one row of Si atoms 
along the [110] direction. Various structures of adsorbed metals and their electronic properties are 
examined. It is found that Al and Ga atomic wire structures with metallic property are strongly 
unstable towards the formation of buckled metal dimers leading to semiconducting behavior. Indium 
atomic wire, however, displays only marginal preference towards the formation of symmetric dimers 
staying close to the metallic limit. The reasons behind the lack of metallic atomic wires are explored. 
In addition, a direction is proposed for the realization of metallic wires on the dihydrogeneted Si(001) . 

PACS numbers: 73.20.-r, 73.21.Hb, 73.90.+f 
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I. INTRODUCTION 

The study of metals on semiconductors dates back to 
the nineteenth century and has seen a vigorous recent re- 
vival due to tremendous interest in Nanotechnology. The 
scanning tunneling microscopy has enabled us to manipu- 
late atoms, place them at will on different surface sites to 
create exotic artificial atomic scale structures with novel 
electronic propertiesiii2i2ii The placement of metal atoms 
such as Al, Ga and In on Si(001) may lead to the forma- 
tion of low-dimensional structures^, exhibiting significant 
new electronic and transport properties. Atomic scale 
structures themselves have technological applications in 
developing atomic scale devices*^ In particular, realiza- 
tion of a one dimensional metallic nanowire is of great 
importance because of its possible use as metallic inter- 
connect in nano-devices. 

There is much current activity in bottoms up approach 
where free standing atomic and nanowires for a variety 
of atoms, e.g., K, Al, Cu, Ni, Au and Si have been 
studiediiSiiiiSiiiiiSii&ii. The geometrical structures of 
such free standing wires and their electronic properties 
have been discussed. A general finding is that a zigzag 
structure in the form of an equilateral triangle is most 
stable£>SiUil4. This can be understood as arising primar- 
ily due to the maximization of coordination number for 
each atom in a quasi ID structure. Another structure, a 
local minimum on energy surface, but not terribly stable, 
is a wide angle isosceles triangle which somehow is rem- 
iniscent of the bulk environment. For example, Si which 
is four fold coordinated in the bulk (tetrahedral angle 
~ 109°) showsi^ a local minimum at an angle of ~ 117°. 
In general, free standing atomic wires tend to be metallic 
(have bands crossing the Fermi level) but these wires in 
practice have to be supported. Silicon is the most widely 
used substrate for practical applications and the low in- 
dex surfaces, Si(001) is the surface of choice. With the 
downward spiral toward nano devices, it is desirable to 
investigate the electronic properties at the lowest pos- 
sible coverages. It is in this context that the study of 



metals like Al, Ga and In at submonolayer coverages on 
Si(001) take on the added importance. The interaction 
of metal nanowires with substrate can significantly alter 
the electronic properties, and not always in the desired 
direction. 

There continues to be a persistent search for metal- 
lic nanowires on clean Si(001): 2x 1 surface as 
well as on the hydrogen terminated Si(001): 2x1 
surfaco 4 i 15 i 16 i 17 i 18 i 19 i 20 . In an early studjsi^, Batra pro- 
posed the formation of a zigzag atomic chain of Al on 
Si(001): 2x1 but energetically it was not the most fa- 
vored structure. In recent studies, it was shown that the 
zigzag Al chain is hard to fabricate as it is energetically 
1.6 eV higher than the most favorable structure. How- 
ever, this chain does not undergo a Peierls distortion and 
remains semimetallic in characten^. 

Recently, the hydrogen terminated Si(001) has become 
one of the surfaces of choice for growing atomic scale 
structures. The hydrogen terminated Si(001) can have 
various reconstructed patterns such as 2 x 1, 3 x 1 and 
lxl depending on the hydrogen coverage and the exper- 
imental environment! 21 ' 22 : 2 ?' i 2 ^i 2 ^i 2 7 Watanabe et alm& 
explored the growth of Ga on the patterned monohy- 
dride Si(001): 2x 1. Using STM tip a monohydried 
Si(001) may be patterned by removing hydrogen atoms 
on a chosen row of surface Si atoms either along [110] 
or [110] direction. Watanabe et a L 19 ^ 20 examined several 
possible structures of Ga on such a patterned monohy- 
dried Si(001): 2x1 surface and they found one dimen- 
sional structures made of small Ga clusters. However, 
the structures turned out to be either semimetallic or 
semiconducting. 

Recently it has been shown in an experiment that an 
ideal hydrogen terminated Si(001): lxl surface can be 
achieved by wet-chemical etching^ Furthermore, it is 
known that with the help of a STM tip, some selected 
hydrogen atoms from the surface may be desorbed. We 
therefore, explore the adsorption of Al, Ga and In on a 
patterned dihydrogeneted Si(001): lxl where Si atoms 
on a single row along the [110] direction are depassivated. 
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In other words, the Si atoms on a single row along the 
[110] direction have two dangling bonds each while all 
other Si atoms on the surface are saturated with hydro- 
gen atoms. Various possible structures for Al, Ga and In 
on the dihydrogeneted Si(001): lxl and their properties 
are studied and a direction is proposed for the realization 
of a metallic atomic wire on the dihydrogeneted Si(001): 
lxl. 




FIG. 1: Silicon and hydrogen atoms within the supercell 
(3x2) for the patterned dihydrogeneted Si(001):lxl. The 
large and small circles represent Si and hydrogen atoms re- 
spectively. The top layer Si atoms on the rightmost row have 
two dangling bonds each as they are not passivated by hydro- 
gen atoms. 



II. METHOD 

First principle total energy calculations were carried 
out within the density functional theory at zero tem- 
perature using the VASP code2£. The wave functions 
are expressed by plane waves with the cutoff energy 
\k + G\ 2 < 250 eV. The Brillouin Zone (BZ) integrations 
are performed by using the Monkhorst-Pack scheme with 
4x4x1 fc-point meshes for 3x2 primitive cells. Ions are 
represented by ultra-soft Vanderbilt type pseudopoten- 
tials and results for fully relaxed atomic structures are 
obtained using the generalized gradient approximation 
(GGA). The preconditioned conjugate gradient method 
is used for the wave function optimization and the con- 
jugate gradient method for ionic relaxation. 



BAND STRUCTURE DOS 




r x m r y 10 20 30 



FIG. 2: Band structure (left panel) and density of states 
(right panel) for the patterned dihydrogeneted Si(001) sur- 
face shown in Fig. The dotted line represents the fermi 
level. The V, X, M and Y points are shown in the inset of the 
right panel. 

The Si(OOl) surface is represented by a repeated slab 
geometry. Each slab contains five Si atomic planes. The 
bottom layer Si atoms are passivated by hydrogen atoms. 
In addition, within the supercell, two consecutive rows of 
Si atoms extending along the [110] direction on the top 
layer are passivated with Hydrogen atoms (see Fig. 
Consecutive slabs are separated by a vacuum space of 9 
A. The Si atoms on the top four layer of the slab and 
hydrogen atoms attached to top layer Si atoms are al- 
lowed to relax while Si atoms in the bottom layer of the 
slab and the passivating Hydrogen atoms are kept fixed 
to simulate the bulk like termination. The convergence 
with respect to the energy cutoff and the number of k 
points for similar structures has been examined earliest. 



III. RESULTS AND DISCUSSIONS 

The lowest energy structure of the patterned dihydro- 
geneted Si(001) surface is shown in Fig. ^ The surface 
retains its unit periodicity along the y[110] direction. 
The unsaturated Si atoms on the surface (third row of 
Si atoms on the top layer in Fig. ^) forms a dangling 
bond wire extending along the y direction. Figure El 
shows the band structure (left panel) and the density of 
states (right panel) for the dangling bond wire. A wide 
band gap (~1.3 eV) around the fermi level indicates that 
the surface is semiconducting in nature. The band gap is 
reflected in the density of states plot with vanishing den- 
sity of states around the fermi level. This is in contrast 
to the metallic nature of dangling bond wire on monohy- 
dride Si(001): 2x1 surfacei 19 ' 20 The reason behind the 
non-metallic nature of the dangling bond wire on the di- 
hydrogeneted Si(001) is that the unit cell has two free 
electrons and they are fully accommodated in a single 
band below the fermi level. 
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Metal atoms are adsorbed on the surface shown in 
Fig. ^ to examine the possibility of the formation of a 
metallic nanowire supported on the Si substrate. Exper- 
iments have shown that the Al and Ga atoms can diffuse 
easily^*! on the hydrogen terminated surface and there- 
fore, the exposed Al, Ga and In atoms are expected to 
diffuse and nucleate around the Si dangling bonds on 
the surface. Thus it may be possible to form nanowire 
of metal atoms supported on the substrate. We are in 
search of a supported atomic wire that will be metallic 
in character. The adsorption of metals like Al, Ga and 
In is studied as a part of this search. 
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FIG. 3: Schematic representation of the atomic positions on 
the top of the slab within the supercell. Small circles represent 
the hydrogen atoms on the top of the slab and large circles 
represent the fiest layer Si atoms. The Si atoms with dangling 
bonds are marked as 1, 2 and 3 respectively. Fhe sites Si and 
S2 have identical surrounding and similarly sites Pi and P2 
have identical environment. 
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FIG. 4: The band structure (left panel) and the density of 
states (right panel) corresponding to the S1S2 configuration 
of Al. The dotted line shows the fermi level. 

Here we consider the Al adsorption on the patterned 
dihydrogeneted Si(OOl) surface shown in Fig. ^ The 
surface offers various possible sites for the Al adsorp- 
tion. Two different kinds of sites are denoted by S (Si 
or S2) and P (Pi or P2) respectively as shown in Fig. [3] 




FIG. 5: The band structure corresponding to the P1P2 con- 
figuration of Al. The fermi level is indicated by the dotted 
line. 






FIG. 6: This corresponds to the most favorable for the Al- 
dimer configuration: (A) The charge density plot on the y- 
z plane including the Al atoms and (B) the band structure 
where the dotted line indicates the fermi level. 



The sites vertically above the hydrogen free Si atoms (Si 
atoms with dangling bonds) are denoted as T (Ti or 
T2). The configurations considered here are T1T2 (one 
Al is placed on top Si atom marked as 1 and the other 
Al atom is placed on top of the Si atom marked as 2), 
S1S2 (one Al is placed at Si site and the other one is 
placed at S2 site) and P1P2 (one Al is placed at Pi site 
and the other is placed at P2 site) respectively. For the 
T1T2 configuration, the Al atoms are allowed to move 
along the z direction only, this configuration is found to 
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be least favorable. The total energies of other configura- 
tions are calculated with respect to the total energy for 
the TiT 2 configuration. The S1S2 and P1P2 configura- 
tions are more favorable than the T1T2 by 0.42 and 0.71 
eV respectively. The S1S2 configuration is interesting be- 
cause the band structure and the density of states (see 
Fig- El f° r this structure indicate the metallic behavior of 
Al atomic wire extending along the y direction. On the 
other hand, the Al structure with the P1P2 configuration 
shows non-metallic behavior (see a band gap around the 
fermi level in Fig. 0. For the S1S2 configuration, each 
Al atom forms two bonds with two Si atoms (each Si-Al 
bond length is ~ 2.6 A ) and the surface retains its unit 
periodicity along the y direction. Therefore, the single 
free electron of the Al atom in the unit cell is respon- 
sible for the partially filled bands crossings through the 
fermi level. Consequently the density of states increases 
around the fermi level and the atomic wire (with S1S2 
configuration) becomes metallic in character. For the 
P1P2 configuration, we note that Al atoms make strong 
bonds with Si atoms (bond length ~ 2.5 A while they fail 
to make a bond among themselves (distance between two 
Al atoms ~ 3.2 A). In this configuration, the periodic- 
ity along y direction is doubled compared to that for the 
S1S2 configuration. Two electrons from two Al atoms be- 
come non-itinerant in the unit cell, prefer to be occupied 
by a single band and hence the Al structure (with P1P2 
configuration) extending along the y direction becomes 
non-metallic. Similar to the zigzag Al chain considered 
earlie»i5ii£ on the bare Si(001): 2x1, we also considered a 
zigzag Al chain configuration by displacing the Al atom 
at Si site by A along the +ve x direction and the Al 
atom at S2 site by A along the -ve x direction. This 
zigzag chain configuration turned out to be unstable and 
readily reverted to the S1S2 configuration. 
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FIG. 7: The band structure (left panel) and the density of 
states (right panel) corresponding to the S1S2 configuration 
of Ga. The fermi level is shown by the dotted line. 

However, the most favorable configuration is the one 
where the Al atom at Si site is shifted along the positive y 
direction by 0.2 A while that at the S2 is shifted along the 
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FIG. 8: This corresponds to the most favorable configuration 
of Ga (Ga-dimer configuration): (A) The charge density plot 
on the y-z plane including the Ga atoms and (B) the band 
structure where the fermi level is indicated by the dotted line. 

negative y direction by 1.6 A to form a buckled Al dimcr 
with a bond length of ~ 2.6 A. The buckled Al dimers 
can be seen in charge density plot in the Fig.EjA) where 
the large circles represent Al atoms while the small circles 
represent Si atoms (just below the Al wire). We will call 
this configuration as Al-dimer configuration and has a 
behavior akin to buckled Si dimer on Si (001). This is 
clear from the charge density plot in Fig.^A). The total 
energy of this configuration is -0.94 eV compared to the 
T1T2 configuration, i.e., this Al-dimer configuration is 
favorable over S1S2 configuration by 0.52 eV. This energy 
is gained mostly due to the optimization of Al-Si bonds 
and the superiority of Al-Al bond with respect to Al-Si 
bond. The nature of the Al wire extending along the y 
direction becomes non-metallic as can be seen from the 
band structure with an energy gap around the fermi level 
(see Fig. E{B)). We therefore, conclude that the metallic 
Al nanowire ( corresponding to the configuration S1S2 
on the dihydrogeneted Si(001) can not be achieved in 
practice. 

We nest consider the Ga adsorption on the same pat- 
terned dihydride Si(001). It is known that at some cov- 
erages, Ga behaves in a different way compared to Ali&. 
The T1T2 configuration is the least favorable one. How- 
ever, unlike the case of Al, the S1S2 configuration for Ga 
(total energy -0.41 eV compared to T X T 2 configuration) 
is slightly more favorable by 0.02 eV compared to the 
P1P2 configuration. The Ga atomic wire with the S1S2 
configuration is metallic with a peak for density of states 
around the fermi level (see Fig. 0). However, the most 
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stable configuration is very similar to that we found for 
Al. Two Ga atoms form buckled dimer as can be seen 
from the charge density plot in Fig.|SjA) (the large circles 
represent the Ga atoms and the small circles represent Si 
atoms just below the Ga atoms). This most favorable 
configuration is denoted as Ga-dimer configuration. The 
atomic wire made of buckled Ga dimers extending along 
the y direction is non-metallic in nature and this can be 
seen from the band structure plot in Fig.[SlB). Similar to 
the Al case, we therefore, conclude that stable metallic 
Ga wire can not be realized on patterned dihydrogenetd 
Si(OOl). We however, notice that the total energy differ- 
ence between the Ga-dimer configuration and the S1S2 
configuration is 0.41 eV which is lower by 0.1 eV com- 
pared to that for the Al case. The reduction in the total 
energy difference between the S1S2 and dimer configura- 
tion for Ga encourages us to study the In adsorption on 
the same patterned dihydrogented Si(001). 

Unlike the case of Al and Ga the P1P2 is the least favor- 
able configuration for In. The total energy for this config- 
uration compared to the T1T2 configuration is +0.08 eV. 
The total energy for the S1S2 configuration compared to 
the T1T2 configuration is -0.20 eV. The In atomic wire 
with the S1S2 configuration is metallic in character with 
large density of states around the fermi level and this can 
be seen from the band structure (left panel) and density 
of states (right panel) plot in Fig. EI The most favorable 
configuration for In is a non-buckled In dimer as seen 
from the charge density plot in Fig. HOf Ai where the In 
atoms (large circles in the figure) at Si and S2 positions 
move towards each other by 0.3 A to form a weak dimer. 
Consequently, the nanowire consisting of unbuckled weak 
In dimers becomes non-metallic (see the band structure 
in the Fig. ^3 (B)). Wc note that the total energy for 
this configuration is -0.27 eV which is more favorable to 
the S1S2 configuration only by 0.07 eV. At room tem- 
perature the thermal energy is expected to be sufficient 
to break these weak dimer bonds leading to metallic be- 
havior of In wires on the dihydrogenetd Si(001). Here we 
add that there are situations where the experimentally 
observed structure corresponds to some local minimum 
energy structure^ffi Therefore, an experiment on this 
system is desirable to confirn if In atomic wire can be 
realized. 



IV. CONCLUSION 

First principle electronic structure calculations are per- 
formed to examine the possibility for the formation of sta- 
ble metallic atomic wires on the dihydrogenetd Si(001). 
Adsorption of metals like Al, Ga and In are considered for 
this purpose. We found that the Al and Ga nanowire con- 
figurations with metallic character are strongly unstable 
towards the formation of buckled metal dimers leading 
to semiconducting behavior. However, the metallic In 
wire is weakly unstable because the total energy corre- 
sponding to the metallic wire configuration is very close 
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FIG. 9: The band structure (left panel) and the density of 
states (right panel) corresponding to the S1S2 configuration 
of In. The fermi level is represented by the dotted line. 
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FIG. 10: This corresponds to the most favorable configura- 
tion for In: (A) The charge density plot on the y-z plane 
including the In atoms and (B)The band structure, dotted 
line represents the fermi level. 



to the most favorable non-metallic (weakly dimerised) 
configuration. The thermal energy may be able to break 
the weak bonds between In atoms and thus a metallic 
In wire may be realized on Si(001). Our results clearly 
indicates that as we go from Al to In via Ga, the metal- 
lic nanowire configuration approaches towards the most 
favorable one. We are hopeful that this work will encour- 
age further experimental studies of In atomic wires on a 
patterned dihydrogeneted Si(001). 
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